Pairing fluctuation theory of high T c superconductivity in the presence of nonmagnetic impurities 
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The pseudogap phenomena in the cuprate superconductors requires a theory beyond the mean field BCS level. 
A natural candidate is to include strong pairing fluctuations, and treat the two-particle and single particle Green's 
functions self-consistently. At the same time, impurities are present in even the cleanest samples of the cuprates. 
Some impurity effects can help reveal whether the pseudogap has a superconducting origin and thus test various 
theories. Here we extend the pairing fluctuation theory for a clean system [Phys. Rev. Lett. 81, 4708 (1998)] 
to the case with nonmagnetic impurities. Both the pairing and the impurity T matrices are included and treated 
self-consistently. We obtain a set of three equations for the chemical potential fi, T c , the excitation gap A(T C ) 
at T c , or /i, the order parameter A sc , and the pseudogap A pg at temperature T < T c , and study the effects 
of impurity scattering on the density of states, T c and the order parameter, and the pseudogap. Both T c and 
the order parameter as well as the total excitation gap are suppressed, whereas the pseudogap is not for given 
T < T c . Born scatterers are about twice as effective as unitary scatterers in suppressing T c and the gap. In 
the strong pseudogap regime, pair excitations contribute a new T 3 ^ 2 term to the low T superfluid density. The 
initial rapid drop of the zero T superfluid density in the unitary limit as a function of impurity concentration m 
also agrees with experiment. 



PACS numbers: 74.20.-z, 74.25.Fy, 74.25.-q 



I. INTRODUCTION 



The pseudogap phenomena in high T c superconductors 
have been a great challenge to condensed matter physicists 
since over a decade ago. These phenomena manifestly con- 
tradict BCS theory by, e.g, presenting a pseudo excitation gap 
in single particle excitation spectrum. Yet the origin of the 
pseudogap and, in general, the mechanism of the supercon- 
ductivity are still not clear. Many theories have been pro- 
posed, which fall into two classes, based on whether the pseu- 
dogap has a superconducting origin. Some authors propose 
that the pseudogap may not be related to the superconductiv- 
ity; instead, it is associated with another ordered state, such 
as the antiferromagnetism related resonating valence bond 
(RVB) state,!!! d-density waveB and spin density wave orderB 
On the other hand, many others believe that the pseudogap 
has the same origin as the superconductivity, such as the 
phase fluctuation scenario of Emery and Kiv^knriLand the 
various precursor superconductivity scenarios. EEoffilo Previ- 
ously, Chen and coworkers have worked out, within the.ru 
cursor conductivity school, a pairing fluctuation theoryEniil 
which enables one to calculate quantitatively physical quanti- 
ties such as the phase diagram, the superfluid density, etc. for 
a clean system. In this theory, two-particle and one-particle 
Green's functions are treated on an equal footing, and equa- 
tions are solved self-consistently. Finite center-of-mass mo- 
mentum pair excitations become important as the pairing in- 
teraction becomes strong, and lead to a pseudogap in the exci- 
tation spectrum. In this context, these authors have been able 
to obtain a phase diagram and calculate the superfluid density, 
in (semi)quantitative agreement with experiment. 

However, to fully apply this theory to the cuprates, we need 
to extend it to impurity cases, since impurities are present even 
in the cleanest samples of the high T c materials, such as the 
optimally doped YBa2Cu307_s (YBCO) single crystals. In 
addition, this is necessary in order to understand the finite fre- 
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quency conductivity issue. Furthermore, study of how various 
physical quantities respond to impurity scattering may help 
to reveal the underlying mechanism of the superconductivity. 
For example, it can be used to determine whether the pseudo- 
gap has a superconducting originJlj Particularly, it is impor- 
tant to address how T c and the pseudogap itself vary with im- 
purity scattering, especially in the underdoped regime. To this 
end, one needs to go beyond BCS theory and include the pseu- 
dogap as an intrinsic part of the theory. Due to the complexity 
and technical difficulties of this problem, there has been vir- 
tually no work in the field on this important problem. 

Among all physical quantities, the density of states (DOS) 
N(ui) close to the Fermi level (u> = 0) is probably most sensi- 
tive to the impurities. Yet different authors have yielded con- 
tradictory results in this regard. BCS-based imruirity T ma- 
trix calculations predict a finite DOS at uj = OjliO which has 
been used to explain the crossover from TLto T 2 power law 
for the low temperature superfluid densityJij Nonperturbative 
approaches have also been_studied and have yielded different 
results. Senthil and FishenlJ find that DOS,vanishes according 
to universal power laws, Pepin and LeefJ predict that N(lj) 
diverges as w ^ 0, assuming a strict particle-hole symme- 
try, and ZiegleiU and coworkers' calculation shows^a rigorous 
lower bound on N(u>). Recently, Atkinson et a/.Ej try to re- 
solve these contradictions by fine-tuning the details of the dis- 
orders. Nevertheless, all these calculations are based on BCS 
theory and cannot include the pseudogap in a self-consistent 
fashion, and thus can only be applied to the low T limit in 
the underdoped cuprates. Therefore, it is necessary to extend 
the BCS-based calculations on impurity issues to include the 
pseudogap self-consistently. 

In this paper, we extend the pairing fluctuation theory from 
clean to impurity cases. Both the impurity scattering and 
particle-particle scattering T matrices are incorporated and 
treated self-consistently. This goes far beyond the usual self- 
consistent (impurity) T-matrix calculations at the BCS level 
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by, e.g., Hirschfeld and others. In this context, we study 
the evolution of T c and various gap parameters as a function 
of the coupling strength, the impurity concentration, the hole 
doping concentration and the impurity scattering strength. In 
addition, we study not only the Born and the unitary limits, but 
also at intermediate scattering strength. We find that the real 
part of the frequency renormalization can never be set to zero, 
the chemical potential adjusts itself with the impurity level. 
As a consequence, the positive and negative strong scattering 
limits do not meet. The residue density of states at the Fermi 
level is generally finite at finite impurity concentrations, in 
agreement with what has been observed experimentally. Both 
T c and the total excitation gap decrease with increasing impu- 
rity level, as one may naively expect. Born scatterers are about 
twice as effective as unitary scatters in suppressing T c and the 
gap. In the unitary limit, the zero temperature superfluid den- 
sity decreases faster with rii when n, is still small, whereas in 
the Born limit, it is the opposite. At given T < T c , both the 
order parameter and the total gap are suppressed, but the pseu- 
dogap is not. Finally, incoherent pair excitations contribute an 
additional T 3 / 2 term to the low T temperature dependence of 
the superfluid density, robust against impurity scattering. 

In the next section, we first review the theory in a clean sys- 
tem, and then present a theory at the Abrikosov-Gor'kov level. 
Finally, we generalize it to include the full impurity T-matrix, 
in addition to the particle-particle scattering T-matrix, in the 
treatment, and obtain a set of three equations to solve for \x, T c 
and various gaps. In Sec. we present numerical solutions 
to these equations. We first study the effects of impurity scat- 
tering on the density of states, then study the effects on T c and 
the pseudogap at T c , followed by calculations of the effects 
on the gaps and the superfluid density below T c . Finally, we 
discuss some related issues, and conclude our paper. 



II. THEORETICAL FORMALISM 



A. Review of the theory in a clean system 

The cuprates can be modeled as a system of fermions which 
have an anisotropic lattice dispersion = 2tn(2 — cosk x — 
cos k y ) + 2t±(l — cos k z ) — //, with an effective, short range 
pairing interaction Vk,k' = .g^kVks where g < 0. Here tn 
and t± are the in-plane and out-of-plane hopping integrals, 
respectively, and [i is the fermionic chemical potential. For 
the cuprates, t± <C £11. The Hamiltonian is given by 

W° = X^^i^ 

kcr 

+ X t/ k,k'C,j < . +q/2T c^ k+q/2i c_ k , +q/2i c k , +q/2t (,l) 

kk'q 

The pairing symmetry is given by ip^ = 1 and (cosfc^ — 
cos k y ) for s- and d-wave, respectively. Here and in what fol- 
lows, we use the superscript "0" for quantities in the clean 
system, to be consistent with the notations for the impu- 
rity dressed counterpart below. For brevity, we use a four- 
momentum notation: K = (k, iui), J2k ~ T J2k u> etc - 

To focus on the superconductivity, we consider only the 
pairing -channel, following early work by Kadanoff and 
Martin;Ej the self-energy is given by multiple particle-particle 
scattering. The infinite series of the equations of motion are 
truncated at the three-particle level G3, and G3 is then factor- 
ized into single- (G) and two-particle (G2) Green's functions. 
The final result is given by the Dyson's equations for the sin- 
gle particle propagator [Refer to Refs. (|9}[l2|) for details] 



X°(K) = G° 1 {K)-G° \K) 



Q 



t\Q)Gl{Q-K)yl_ 



q/2 



and the T matrix (or pair propagator) 

t°(Q) = t° sc (Q)+t° pg (Q) 

with 



(2a) 



(2b) 



The excitation gap forms as a consequence of Cooper pair- 
ing in BCS theory, while the superconductivity requires the 
formation of the zero-momentum Cooper pair condensate. As 
these two occurs at the same temperature in BCS theory, one 
natural way to extend BCS theory is to allow pair formation 
at a higher temperature (T*) and the Bose condensation of the 
pairs at a lower temperature (T c ). Therefore, these pairs are 
phase incoherent at T > T c , leading to a pseudogap without 
superconductivity. This can nicely explain the existence of 
the pseudogap in the cuprate superconductors. Precursor su- 
perconductivity scenarios, e.g., the present theory, provides a 
natural extension of this kind. At weak coupling, the contri- 
bution of incoherent pairs is negligible and one thus recov- 
ers BCS theory, with T* — T c . As the coupling strength 
increases, incoherent pair excitations become progressively 
more important, and T* can be much higher than T c , as found 
in the underdoped cuprates. In general, both fermionic Bo- 
goliubov quasiparticles and bosonic pair excitations coexist at 
finite T <T C . 



t° sc (Q) = —^S(Q) , 



T 



where A sc = at T > T c , and 

C(Q) 



l + 9X°(Q) 



(2c) 



(2d) 



where A sc is the superconducting order parameter, Gq(K) = 
l/(iu> — £k) is the bare propagator, and 



X °(Q) = X G°(K) G° (Q - K) p k _ q/2 (2e) 

K 



is the pair susceptibility. This result can be represented dia- 
grammatically by Fig. [l|. The single and double lines denote 
the bare and full Green's functions, respectively, and the wig- 
gly double lines denote the pair propagator. 

The superconducting instability is given by the Thouless 
criterion 



1 + .9X» = 0, (T<T C ), 



(3) 
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FIG. 1: Diagrams for the Dyson's equations in a clean system, 
which leads to the approximation 



where the pseudogap is defined by 



A; 



E *™«) 



(5) 



As a consequence, the self-energy takes the standard BCS 
form 



(6) 



where A 2 = Aj c + A 2 pg , and A k = Atp k . In this 
way, the full Green's function G°(K) also takes the stan- 
dard BCS form, with the quasiparticle dispersion given by 
E k = \J e 2 . + A 2 ^ 2 .. So does the excitation gap equation 



i+sE 



l-2/f.Ek) 9 



(7) 



X' 



(a) 



vV 

(b) 



(c) 



FIG. 2: Examples of (a) bridging, (b) crossing diagrams, and (c) 
higher order terms neglected in the calculation. 



where in the real space the impurity term is given by 

Hi = J dx u(x - x i )i/> t (x)?/>(x) , 



(10) 



with k(x) = u5(x) for isotropic s-wave scattering. 

To address the impurity scattering, we begin at the 
Abrikosov-Gor'kov (AG) level,E3E3£j which is a good ap- 
proximation in the Born limit. Following AG, we include 
all possible configurations of impurity dressing, but exclud- 
ing bridging diagrams like Fig. ||(a), crossing diagrams like 
Fig. ^(b), and higher order terms like Fig. |](c). The dashed 
lines denote impurity scattering, and the crosses denote the 
impurity vertices. We clarity, in most diagrams, we do not 
draw the fermion propagation arrows. It is understood, how- 
ever, that they change direction at and only at a pairing vertex. 
As in subsec. II A, we use plain double lines to denote the 
Green's function (G°) fully dressed by the pairing interaction 
but without impurity scattering, i.e., 



G°{K) = 



1 



where 



q/2 



(ID 



(12) 



We emphasize that although this equation is formally identi- 
cal to its BCS counterpart, the A here can no longer be in- 
terpretted as the order parameter as A pg ^ in general. For 
self-consistency, we have the fermion number constraint 



= 2J2G°(K)=2J2 



K 



(8) 



The gap equation Eq. (0), the number equation Eq. ([$]), and 
the pseudogap parametrization Eq. (|j) form a complete set, 
and can be used to solve self-consistently for T c , fi(T c ), and 
A(T C ) by setting A sc = 0, or A(T), and A pg (T) at 

given T < T c . Here v 2 = \ (1 - e k /E k ), as in BCS. 



B. Impurity scattering at the Abrikosov-Gor'kov level 

For simplicity, we restrict ourselves to nonmagnetic, elas- 
tic, isotropic s-wave scattering. At the same time, we will 
keep the derivation as general as possible. In the presence of 
impurities of concentration nu the Hamiltonian is given by 



H = H° + H! 



(9) 



However, since we will address the impurity dressing of the 
pairing vertex or, equivalently, the pair susceptibility x(Q), 
we assume the pair propagator t(Q) in the above equation is 
already dressed with impurity scattering, with 



t(Q)=t sc (Q) + t pg {Q) , 



and 



A 



tsc{Q) 



E wq) 



T 



-s(Q) , 



9X(Q) ' 



(13) 



(14) 



(15) 



as in the clean case. The shaded double lines denote impurity- 
dressed full Green's function G, and "shaded" single lines de- 
note impurity-dressed bare Green's function (which we call 
Go), i.e, 



G (if) 



iui - e k - G Q , 



(16) 
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r (K\Q) 
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FIG. 3: Feynman diagrams for the impurity dressed full Green's function. 



where the bar denotes impurity average Gq^ = 
n.j |u(k — k')\ 2 Go(K'). We use open circles to de- 
note bare pairing vertex "f(K\Q) = <Pk+q/2, and shaded 
circles full pairing vertex T(K\Q), where Q is the pair 
four-momentum. To obtain the Feynman diagrams for 
the impurity dressed full G, we first expand the pairing 
self-energy diagram as an infinite series which contains only 
bare single particle Green's function and pair propagators, 
and then insert all possible impurity scattering on the single 
particle propagators at the AG level. We assume that the pair 
propagators are always self-consistently dressed by the im- 
purity scattering. After regrouping all non-impurity dressed 
lines on the left, the final result for the diagrams is shown 
in Fig. [3](a). Here following AG, the subdiagrams inside 
the two impurity legs are assumed to be self-consistently 
dressed by impurity scattering. To make direct comparison 
with the BCS case easier, we present the corresponding 
diagrams for the BCS case in Appendix |X| The first term on 
the right hand side (RHS) of Fig. ||(a) contains all diagrams 
without impurity dressing (except via the pair propagators). 
The second term corresponds to the third term of the first 
equation in Fig. The third term corresponds to the last 
term, the last term to the second. The fourth and the fifth 
together correspond to the fourth term [see Fig. |(a)]. The 
fifth term in Fig. |](a) arises since the two impurity legs can 



cross two separate pairing self-energy dressing parts; it can 
be eliminated using the equality shown in Fig. ||(b). Here 
the shaded elliptical region denotes self-consistent impurity 
dressing of the double pairing vertex structure inside the 
two impurity scattering legs, as shown in Fig. |](a). It is 
worth pointing out that these diagrams reduce to their BCS 
counterpart if one removes the pairing propagators. The 

Dyson's equation Gq 1 ^) — G[j 1 (K) — Gq jUJ can then be 
used to eliminate the fourth term in Fig. ||(a). We now obtain 
the greatly simplified diagrams for G as shown in Fig. [|(c), 
which can be further reduced into Fig. ||(e), upon defining a 
reducedpairing vertex Tq(K\Q) (shaded triangles) as shown 
inFig.gd), 



T (K\Q) = 

+n z Kk - k')| 2 G (Q ~ K')T(K'\Q)G(K') . (17) 



One can then read off the impurity-induced quasiparticle self- 
energy E'(K) = G° (if)— G -1 (iT) immediately, as shown 
inFig.j^e), 



(18) 
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FIG. 4: Feynman diagrams for pairing vertex. 



where the impurity average 

G„ = n,^Kk-k')| 2 G(F') , 

k' 

and the "full" self-energy 



(19) 



E(if ) = ^ r (F|Q)t(Q)r(F|Q)G (Q - if) . (20) 
Q 

Therefore, we have finally 

G-\K) = G ~\K) -E'(if) ^iuj-e^-Gu - E(if ) 
= »w - e k - E(if) = G^F) - E(if) , (21) 

where we have defined the renormalized frequency icl> = iuj — 
G w and the "bare" Green's function 



G (K) 



1 



IUJ — £ k 



(22) 



Note here G (K) + G Q (K). 

Now we deal further with the pairing vertex T(K\Q). First, 
we notice that the impurity-dressed double-vertex structure in 
Fig. 0(a) can be simplified as shown in Fig. 0(a) using the 
approximation 



53*(Q)/(jr,Q)w £t(Q) 



/(^,0) , 



(23) 



where /(if, Q) is an arbitrary slow-varying function of Q. 
This is due to the fact that t(Q) diverges as Q — > at T < 
T c . The Dyson's equation for the Green's function G°(— F) 
in a clean system is also used in getting the second line of 
Fig. 0(a). Therefore, we have approximately the impurity- 
dressed pairing vertex as shown in Fig. [|(b), which implies 
the equality shown in Fig. 0(c). Fig. 0(c) can be written as 



G (K )T(K \Q)G(Q - K ) = 



r (K\Q)G(Q ~ K) 
Gq (if) — G w + Grj )U 



= G (K)T (K\Q)G(Q - K) . (24) 
This result demonstrates the following important relationship: 



G (K)T(K\Q) = G (K)T (K\Q), (small Q) . (25) 



Using this relationship, now the self-energy can be simplified 
as follows: 

E(if) = Y,T 2 (K\Q)t(Q)G (Q - K) 
Q 

« -A 2 T 2 (K)G (-K) = -A k G (-if) , (26) 

where T (K) = T(K\Q = 0) and A k = AT (K). 

Finally, the impurity dressing of each rung [i.e., x(Q)] °f 
the particle-particle scattering ladder diagrams is topologi- 
cally identical to the impurity dressing of the pairing vertex 
and the two associated single particle lines. And summing up 
all the ladders gives the pairing T matrix. Therefore, the pair 
susceptibility becomes 



X (Q) = ]Tr(F|Q)G (Q-F)G(F)^ k _ q/2 

K 

« r (K\Q)G (Q - if)G(if)^ k _ q/2 (27) 

K 

And the gap equation is given by 
1 + <?X(Q) = 0=1+ <?^r (if)G(F)Go(-if)^ k . (28) 

K 

This result can be easily verified to be consistent with the 
self-consistency condition. Define formally the generalized 
Gor'kov F function: 



Ft (if) = AG (-K)T(K)G(K) 
Using Eq. (p5[), we have 



(29) 



Ft (if) = Ar (F)G (-F)G(F) = A k G (-F)G(if) . 

(30) 

The formal difference between this Ft and that in BCS is that 
(£> k is now replaced by the renormalized vertex Tq(K). One 
immediately sees that the condition 



(31) 



K' 



is consistent with the gap equation Eq. (|28j). However, it 
should be emphasized that the F function so defined does not 
vanish above T c in the pseudogap regime, different from the 
BCS case. 
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FIG. 5: Replacement scheme from the AG level to self-consistent 
impurity T-matrix treatment. 



C. Impurity scattering beyond the AG level 

In this subsection, we include both the impurity scatter- 
ing T-matrix with the particle-particle scattering T-matrix, 
and, thus, go beyond the AG level. We notice that if one 
replaces the second-order impurity scattering subdiagrams at 
the AG level with the corresponding impurity T-matrices, as 
shown in Fig. g[ the derivation for G(K) goes through for- 
mally without modification. Now we only need to determine 
the impurity T-matrices T u and T A t (as well as their com- 
plex conjugate) in terms of their AG-level counterpart, G w 
and Ft = n . \u(k - k')\ 2 F^(K'), respectively. In other 
words, except that iui and Ak now have different expressions, 
everything else remains the same in terms of iu) and Ak (as 
well as their complex conjugate), just as in the BCS case (see 
Appendix Q). 

The Feynman diagrams for T w and T A t are shown in Fig. [| 
To obtain the second line for T A t , we make use of the approx- 
imation in Fig. |](a) to convert the left part of the second and 
the third term on the first line to the full G. This result is 
direct analogy with its BCS counterpart as shown in Fig. |o| 
One can now write down the equations for T w and T A t with- 
out difficulty. 

r w (k,k') = u (k,k') + ^u(k,k")G(x")T w (k",k') 

k" 

+ ^^ !1 (k,k")Gr)Go(Q-i^") 

k" Q 

x Y(K"\Q)t(Q)T Af (K" -Q,K'\Q) ,(32) 

and 

T At (K ~ Q, K'\Q) = J2 «(k" - q, k - q) 

k" 

x G a (Q - X")G(A"')r(A"'|Q)T w (k",k') 
+ t»(k" - q, k - q)G(Q - A"')T At (K" - Q, K'\Q) . 

k" 

(33) 

Note here T A t does not contain the factor A, unlike its BCS 
counterpart. It has the same dimension as T. Both Y and Yq 
now contain the full impurity T-matrix beyond the AG level, 
and the vertex relation Eq. ( j25| ) remains valid. The new ex- 
pression for r is given by 

Y (K\Q) = ^k- q /2 + th t ai (K - Q, K'\Q) . (34) 



So far, we have kept the derivation for a generic elastic scat- 
tering it(k, k'). For isotropic s-wave scattering, u(k, k') = u. 
In this case, T w and T A t are independent of k and k'. Neglect- 
ing the momentum dependence, we obtain 



T a \(uj - Q,lu) 



■ Ek G (Q - K)Y a (K\Q)G(K) 



and 



T 

(35a) 



Q 



X T A t — IjJ) 



(35b) 



Use has been made of the vertex relation Eq. (£5). Here again, 
we need to make use of the approximation Eq. (23|). Defining 

Gu = E k G ( K ) and 

H = L F ^ K ^ = E *To(K)G (-K)G(K) , (36) 



we obtain 

T u = 

and 



u(l - uG- u ) 



1 - uG u 1 - uG 



(37a) 



T A \(u-fl,u>) = 



u 2 j: w G Q (Q - K)Y (K\Q)G(K) 
(l - u&^j (l - uG-J) + u 2 F u M 



1-uGL 



1 - uGo_ u 
Letting il — > 0, the last equation becomes 



AT At H = 



1 - uG u ) ( 1 - uG^ ) + u 2 F u Fi 



(37b) 



(37c) 



The frequency and gap renormalizations are given by 

iuj = iui — Yi u , i& = —iui — , (38a) 



A k = A k + E At 



(38b) 



where E w = n{T u and E A = riiAT A . Here w = (— u>). The 
expression for x(Q) remains the same as in previous subsec- 
tion. 

Ford-wave, T A t = T A = 0, and A k = Ak- Then 
Eq. (37a) is greatly simplified, 



T, = 



1-uGu 



(39) 
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The full Green's function is given by 

iQ — e k 



G(K) 



{iCb - e k )(iu>- e k ) + A*A k 



(40) 



Due to the approximation Eq. (^3j), we are able to bring the 
final result Eqs. (|39|) and (|o|) into the BCS form. It is easy 
to show that they are equivalent to the more familiar form in 
Nambu formalism, as used in Ref. Ol Define 



It should be emphasized, however, that unless u = or 
u = ±oo, the symmetric part of the impurity T-matrix, Ts„, 
can never be set to zero, even if one could in principle have 
Gsui — 0. This means that Ik will always acquire a non- 
trivial, frequency dependent renormalization, £s w . While this 
renormalization is small for weak coupling BCS superconduc- 
tors, it is expected to be significant for the cuprate supercon- 
ductors. 



i(T w + T^), (41) 



and similarly for S w and G u . Here the subscript "A" and "S" 
denote antisymmetric and symmetric part, respectively. Fur- 
ther define 



iu> A =ius - £Uw, £k — e k + ^Suj 
then we obtain (with A* = A) 

iui A + e K 



G{K) 



(iCj A ) 2 - & K - A 



and 



Tau 



Tsui — 



u 2 G A u 



(1-uGsu) 2 -u 2 G% 
u(l - Gsu.) 



(1 - uG Su 



n 2 G 2 



where 



Gau 
Gsu 



E 

k 

E 



IU) A 



(iCja) 2 



e 2 - A 2 



(i& A ) 2 - l 2 K - Al 



(42) 
(43) 

(44a) 
(44b) 

(45a) 
(45b) 



are the antisymmetric and symmetric parts of G u , respec- 
tively. It is evident that Ta and T$ correspond to To and — T3, 
respectively, in the Nambu formalism in Ref. |l5j (and simi- 
larly for Ga and Gs). 



III. NUMERICAL SOLUTIONS FOR d-WAVE 
SUPERCONDUCTORS 

A. Analytical continuation and equations to solve 

Since there is no explicit pairing vertex renormalization for 
d-wave superconductors, i.e., Tq(K) — ip^ or A k = A k , a 
major part of the numerics is to calculate the frequency renor- 
malization. Everything else will follow straightforwardly. 

Numerical calculations can be done in the real frequencies, 
after proper analytical continuation. Since Ts^ ^ 0, S w and 
are independent of each other. To obtain the frequency 
renormalization one has to solve a set of four equations 
for G u , G-u, So;, and self-consistently for given ui. Be- 
cause iQ ^ —iui, and one needs to analytically continue both 
simultaneously, the analytical continuation must be done care- 
fully. For n > 0, iui n — > uj+ = uj + + iS" , and iu>_ n — > = 
= —n < 0, iuj n i — ► ui 1 } = oj- + iTj'!_ and 
— iS" . Here u>± = ±cj — S' ± , and we 
choose uj > and £± > 0. Then we obtain four equations as 
follows: 



uj- — iT," . For ml 



Lr w>0 



G R 
£ 



w<0 
w>0 



E 

k 

E 



j_ - £k 



o;+ + i£" - e k )(u;_ - - e k ) 
w+ - iE" - e k 



(w_ + iL'L - e k )(w+ - tS'; - e k ) + A 2 . 



riiU 



yR 



1 - uGS 

TIM 



1 - «G« 



= E'_ - »T" 



(46) 
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These equations are solved self-consistently for (Ej_, E±), as 
well as the real and imaginary parts of G± u , as a function of 
uj. No Kramers-Kronig relations are invoked in these numer- 
ical calculations. Note in real numerics, we subtract E^ =0 
from E^ so that E^ =0 = 0. This subtraction is compensated 
by a constant shift in the chemical potential /i. 

Having solved the frequency renormalization E w , one can 
evaluate the pair susceptibility in the gap equation Eq. 



where b(x) is the Bose distribution function. 

Numerical solutions confirm that the coefficients a' and b' 
have very weak T dependence at low T. Therefore, in a three- 
dimensional (3D) system, we have A 2 pg ~ T 3 / 2 at low T. As 
the system dimensionality approaches two, the exponent de- 
creases from 3/2 to 1. However, in most physical systems, 
e.g., the cuprates, this exponent is close to 3/2. The prod- 



uct n r , 



jpApg roughly measures the density of incoherent 



X(0) 



E 

K 



(iu> 



Im 



E 

k 



- e k )(zw - 
du> 



^k) + A k 

[1 - 2/M] <& 



A 



where f(x) is the Fermi distribution function. It is easy to 
check that the gap equation Eq. ( ^8| ) does reduce to its clean 
counterpart Eq. (^) as m — * 0. 

The particle number equation becomes 



k 



^(k, W )/(w) = -2Im 

Z7T 



by 



(48) 

The real and imaginary parts of x(Q) m ' e given respectively 



= Im E r ^{G K (w,k)<(fl- W) q. 



x -SI)- /(«,)] 

+ G R (u,k)G A (Sl-u,q-k) 

x [i_/( w )_/( w _n)]} ( 



pairs. 

The gap equatio n (|28| ) [together with Eq. (p7|)], the fermion 
number equation (|48[), and the pseudogap equation ( |5l| ) form 
a closed set, which will be solved self-consistently for T c , /i, 
' and gaps at and below T c . For given parameters T c , \i, and 
gaps, we can calculate the frequency renormalization, E w , and 
(47)then solve the three equations. A equation solver is then used 
to search for the solution for these parameters. The momen- 
tum sum is carried out using integrals. Very densely populated 
data points of oj are used automatically where E w and/or G u 
change sharply. A smooth, parabolic interpolation scheme is 
used in the integration with respect to u>. The relative error of 
the solutions is less than 1.0 x 10~ 5 , and the equations are sat- 
isfied with a relative error on both sides less than 1.0 x 10~ 7 . 
In this way, our numerical results are much more precise than 
those calculated on a finite size lattice. 

The solutions of these equations can be used to calculate 
the superfluid density n s /m. Without giving much details, 
we state here that the impurity dressing of the current vertex 
for the (short-coherence-length) cuprate superconductors does 
not lead to considerable contributions for s-wave isotropic 
scattering and the long wavelength q — > limit. The expres- 
sion for the in-plane n s /m is given formally by the formula 
for a clean system (before the Matsubara summation is carried 



q/2 : 



(49a) 



and 



X "(Sl + iO+,q) 



E 



^Im G r (lj, k)A {Sl - w, q - k) 



x [/(ft- w )-/( w M_ q/2 



(49b) 

where Aq(u>,\s.) = —2ImG R (ui,k) is the "bare" spectral 
function. 

The pseudogap is evaluated via Eq. (|l5|). To this end, we 
expand the inverse T-matrix to the order of SI and q 2 via a 
(lengthy but straightforward) Taylor expansion 

*0+,q) 



= x(tt + *0+,q)- X (0,O) 

= K> + 



)Sl + b'q 2 + ia'{Sl 2 .(50) 



out), as in Ref. 25 



-i = - + P(0) , 
m m 



(52) 



where the in-plane current-current correlation function 
P(Q)Sjj can be simply derived from Eqs. (31) and 



P»(Q) 



(32) of Ref. E5|. The result is given by 



P(Q) = J2 G ( K ) G ( K -Q) 



K 



x G Q (-K)G (Q - K) 
l ^e k _ q / 2 



(AL - A^ k¥ , k _ q 

^ £ k-q/ 2 N 



Oku 



dkn 



(53) 



Here the imaginary part b" vanishes. The term b'q- should 
be understood as b'^q 2 + b' ± q± for a quasi-two dimensional 
square lattice. We keep the imaginary part up to the St 2 order. 
Substituting Eq. ( |5(j( ) into Eq. (fTs]), we obtain 

dSl K + a'{Sl)Sl 



A 



E 



{a'vSl 



b'q 2 ) 2 + « 



a'{Sl) 2 Sl 



;b(Sl) , 

(51) 



Using spectral representation and after lengthy but straightfor- 
ward derivation, we obtain 



= 4ALE 



dui 
2^ 



Im 



(F^,k)) 2 (V6 k ) 



2 ,„2 



~G A ( W) k)F A (u,,k)Ve k 



/M , (54) 
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FIG. 7: Example of (a) the frequency renormalization — an d (b) 
impurity average of Green's function — for a d-wave supercon- 
ductor. The dashed and solid lines denote the real and imaginary 
parts, respectively. The full impurity band for — is shown in the 
inset. The parameters used are \i = 0.92, t±/tu = 0.01, rii = 0.02, 
u = 1, A = 0.08. The energy unit is 4t||, the half bandwidth. 



where 

= G(K)G (-K 



1 



(mD - e k )(iw - e k ) + A k ' 

(55) 

which differs from F(K) by a factor A k . 

As in the clean system, Eq. (|4j) differs from its BCS coun- 
terpart (n s /m) bcs on ly by the overall prefactor, A 2 C , 



m 



A 2 
A 2 



(-) 



m / bcs 



(56) 



For d-wave superconductors, (n s /m) bcs ~ A — BT or A — 
BT 2 at very low T, depending on whether the system is clean 
or dirty. Bearing in mind that A 2 C /A 2 = 1 - A 2 S /A 2 — 
1 - C'T 3 / 2 , we have 

— - A - BT - CT 3/2 , (clean), 
m 

~ A - BT 2 - CT 3/2 , (dirty), 



(57) 



at very low T. In other words, pair excitations lead to a new 
T 3 / 2 term in the low T superfluid density. 



B. Renormalization of the frequency by impurity scattering 
and the density of states 

Except in the Born limit, for u > 4tii, impurity scattering 
usually introduces a sharp resonance close to the Fermi level 
in the frequency renormalization for a d-wave superconduc- 
tor. In addition, it induces an impurity band outside the main 
particle band. Both the low energy resonance and the high 
energy impurity band arise from the vanishing of the real part 
of the denominator of the impurity T-matrix, Eq. (|3^), while 




FIG. 8: Evolution of (a) the frequency renormalization — Im and 



(b) impurity average of Green's function 

.2 



-lmGi 



with u for a d- 

wave superconductor at fixed 7 = riivf — 0.02. A resonance de- 
velops in — Im as u deviates from the Born limit. The weak scat- 
tering (Born limit) is more effective in filling in the DOS within the 
gap. The parameters used are fi = 0.9, t±/tu = 0.01, A = 0.0945. 



the imaginary part is small. In Fig. ^ we show an example of 
(a) the retarded, impurity induced renormalization of the fre- 
quency, — = uj r — oj and (b) the corresponding impurity 
average of the single-particle Green's function, —G^, which 
is related to the density of states by N(u>) = — 2ImG^. The 
curves for — corresponding to the impurity band are re- 
plotted in the inset, to show the strong renormalization of the 
frequency inside the impurity band. (The curves for — has 
been offset so that — E^ =0 = 0. This offset is compensated 
by a shift in the chemical potential /i.) The van Hove singu- 
larity and its mirror image via particle-hole mixing are clearly 
seen in the density of states, and are also reflected in — E^ 
as the small kinks in Fig. ^(a). The real part — ELis usually 
neglected in most non-self-consistent calculations]!!] It is evi- 
dent, however, that it has a very rich structure, and, in general, 
cannot be set to zero in any self-consistent calculations. This 
conclusion holds even in the presence of exact particle-hole 
symmetry, as can be easily told from Eq. (]3^). 

For u < 0, the low energy resonance in — Im E^ will appear 
on the positive energy side in Fig. 0(a). Regardless of the sign 
of u, the resonance peak will become sharper as m decreases 
and as \u\ increases. For larger the resonant frequency 
will be closer to u> = 0, where — Im G R is small because of 
the d-wave symmetry; A smaller n, further reduces — ImGj. 
Both factors help minimize the imaginary part of the denom- 
inator of Eq. ( p^ ) and, thus, lead to a stronger resonance. It 
should be emphasized that a resonance peak does not show 
up in — ImG^ since the resonance in — ImE^ requires that 
— Im G R be small at the resonance point. 

The location of the impurity band is sensitive to the sign 
and strength of impurity scattering. For u < 0, the impurity 
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FIG. 9: Evolution of the DOS N(lu) for a d-wave superconductor (a) 
with u at Tii — 0.02 and (b) with rii in the unitary limit u = 1000. 
There is a dip at w = for small or small The parameters used 
are // = 0.9, t±/ty = 0.01, A = 0.0945. 
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FIG. 10: Residue DOS N(0) at the Fermi level as a function of (a) 
m for u = 1000, 10, 5, 3, 2, and 1, and of (b) u for m = 0.005, 
0.01, 0.02, 0.03, and 0.05. Also plotted in (a) is N(0) as a function 
of y(= m) in the Bom limit. The parameters used are fi = 0.9, 
t±/t\\ = 0.01, A = 0.0945. 



band on the negative energy (left) side of the plot. As \u\ gets 
smaller, the impurity band merges with the main band; as \u\ 
gets larger, it moves farther away, with a much stronger renor- 
malization of u>. For large \u\, the spectral weight under the 
impurity band in Fig. |7](b) is given by 2rii, and the weight in 
the main band is reduced to 2(1 — ni). This leads to a dra- 
matic chemical potential shift as a function of the impurity 
concentration m (as well as u). For u < 0, the impurity band 
will always be filled, so that increasing n>i pushes the system 
farther away from the particle-hole symmetry. For u > 0, on 
the contrary, the impurity band is empty, and the system be- 
comes more particle-hole symmetrical as rii increases from 0, 
and reaches the particle-hole symmetry (in the main band) at 
rii = 1 — n. This fact implies that Pepin and Lee's assump- 
tion of an exact particle-hole symmetry is not justified so that 
their prediction of a diverging N(uj) as uo — > is unlikely to 
be observed experimentally. It should be mentioned that the 
appearance of the impurity band has not been shown in the lit- 
erature, largely because most authors concentrate on the low 
energy part of the spectra only, and do not solve for the full 
spectrum of the renormalization of u> self-consistently. 

In the Born limit, only the product 7 = riiU 2 is a meaning- 
ful parameter, not n, or u separately. In Fig. |[ we plot (a) 
the frequency renormalization -ImS^ and (b) the impurity 
average of Green's function — Im as a function of w for 
various values of u but with a fixed 7 = 0.02. (Note: when u 
is small, this requires an unphysically large m.) In the Born 
limit, these two quantities are identical up to a constant coef- 
ficient. As u increases, a resonance develops in — Im at 
small oj. -Im and -ImG^ become very different. And 
a impurity band develops gradually (it = 0.2 and u = 0.5), 
until it splits from the main band (u = 1). At fixed 7, the Born 
limit is more effective in filling in the DOS in the mid-range of 
u within the gap and smearing out the coherence quasiparticle 
peaks, whereas the large u limit is more effective in filling in 
the DOS in the vicinity of uj = but keeping the quasiparticle 
peaks largely unchanged. In addition, the main band becomes 
narrower at large u than that in the clean system or the Born 
limit, so that part of the spectral weight has now been trans- 
ferred to the impurity band. Also note that the DOS at to = 
is essentially zero in Fig. @(b) because rii is ver y small when 



u becomes large for the current choice of 7 = 0.02. 

The effects of the scattering strength u and the impurity 
density rii on the DOS are shown in Fig. ^(a) and (b), respec- 
tively. For the effect of u in Fig. ^(a), we choose an interme- 
diate rii = 0.02. And for the effect of n,-, we focus on the 
unitary limit, and choose u = 1000. There is a dip in the DOS 
at oj = in both the small u and small m cases, mimicking 
a fractional power law dependence on oj. At higher rii and 
higher u, the DOS is filled in mainly at small oj. 

Shown in Fig. |l^ are the residue DOS at the Fermi level, 
N(0), as a function of (a) the impurity concentration rii for 
different values of u from the unitary limit u = 1000 through 
u = 10, 5, 3, 2, down to 1, and (b) as a function of the scat- 
tering strength u for m = 0.005, 0.01, 0.02, 0.03, and 0.05. 
Figure [l(]( a) indicates that below certain "critical" value of rii, 
N(0) remains essentially zero. This behavior is also implied 
by the presence of the dip at small rii in Fig. |(b). The "crit- 
ical" value for n, in Fig. |o|(a) is clearly scattering strength 
dependent. The smaller u, the larger this value. A replot (not 
shown) of these curves in terms of log 10 iV(0) as a function of 
1 1 rii reveals that for small m, N(0) vanishes exponentially as 
e -A/m^ wnere ^4 i s a constant. For comparison, we also show 
in Fig. [K](a) the Born limit as a function of y(= m). As one 
may expect, the Born limit is rather different from the rest, 
since it is equivalent to a very small u < < 1 and unphysically 
large rii. A similar activation behavior of N(0) as a function 
of u is seen in Fig. |Io|(b), where the "critical" value for u is 
strongly rii dependent. The asymmetry between positive and 
negative u reflects the particle-hole asymmetry at \i = 0.9. It 
should be noted that it is not realistic to vary u continuously 
in experiment. 

An earlier experiment by Ishida et alM suggests that iV(Q) 

1 /2 

varies as n, . In our calculations, however, N(0) does not 
follow a simple power law as a function of n,. The curve 
for u = 1000 in Fig. [uja) fits perfectly with anf — b, with 
a w 0.175, for 0.002 < n t < 0.05. The damping of the zero 
frequency (not shown), — ImE^ =0 , also fits this functional 
form very well, with a « 0.61. The exponents are different 
for different values of u. Our calculation for the rii depen* 
dence of N(0) is consistent with the result of FehrenbacherEZI 
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FIG. 11: Behavior of T c , \i (left inset), and A (right inset) as a func- 
tion of the impurity scattering strength u at n = 0.85, t± /t\\ = 0.01, 
-g/4&\\ = 0.5 and m = 0.05. 



in that A^O, rii) it is strongly u dependent. However, it does 

1/2 

not seem likely that the simple power law N(0) ~ n i / may 
be obtained in an accurate experimental measurement. Fur- 
ther experiments are needed to double check this relationship 
between N(Q) and m. 

From Figs. [7|-|l0| we conclude that for very small n, and u, 
the zero frequency DOS N(0) is exponentially small. At high 
rii, N(0) is finite in both the Born and the unitary limit. How- 
ever, for certain intermediate values of m and u [e.g., u ~ 1 
in Fig. |(a), and m < 0.002 in Fig. |(b)], N{u) vanishes with 
(very small) oj according to some fractional power uj a where 
a < 1. We see neither the universal power laws for N(u>) 
predicted by Senthil and_Fisher,l!Zl nor_the divergent DOS pre- 
dicted by Pepin and LedlJ and others E5 

When the values of u and m are such that N(ui) ~ uu a 
with a < 1, one may expect to see a fractional low T power 
law in the superfluid density. However, such a power law is 
not robust as it is sensitive to the impurity density for a given 
type of impurity. The situation with a negative u is similar to 
Fig. | 

For given chemical potential /j, and the total excitation gap 
A, the calculation of the frequency renormalization with im- 
purities does not necessarily involve the concept of the pseu- 
dogap. It is essentially the "self-consistent" impurity T- 
matrix treatment by Hirschfeld et al\3 except that we now 
have to solve for both the real and imaginary parts of si- 
multaneously in a self-consistent fashion. Our numerical re- 
sults agree with existing calculations in the literature. 

Finally, we emphasize the difference between the self- 
consisteplimpurity T-matrix treatment of the one-impurity 
problerrlBEj and the current many-impurity averaging. For 
the former case, the impurity T-matrix will be given by 

the impurity av- 



t.e. 



Eq. (37c) but with G u replaced by G°. 
erage of the clean G°(K). As a consequence, the position of 
the poles of T w is independent of the renormalized DOS, and, 
therefore, a resonance peak may exist at low u> in the DOS,Eil 
whereas it cannot in the current many-impurity case. 
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FIG. 12: Evolution of (a) A pg (T c ), (b) T c , and /i (inset), as a func- 
tion of the impurity concentration m for both positive (it = 100) 
and negative (u — —100) scattering strength in the unitary limit, 
the Born limit (u = 1, 7 = rii), and intermediate u = —1. Here 
n = 0.85, t±m = 0.01, and -g/4t\\ = 0.5. 



C. Effects of the impurity scattering on T c and the pseudogap 



In this section, we study the influence of impurity scattering 
on the behavior of T c and the pseudogap A pg as a function of 
the coupling strength as well as the hole doping concentration. 

First, we study the effect of the scattering strength u and 
whether it is repulsive (u > 0) or attractive (u < 0). In 
Fig. 1 1 , we plot T c as a function of u, for a pseudogapped 
rf-wave superconductor with m = 0.05. The corresponding 
curves for fi and A = A pg (T c ) are shown in the upper left 



and upper right insets, respectively. All three quantities, T c , 
fi, and A, vary with u. For either sign of u, both T c and A are 
suppressed by increasing \u\. It should be emphasized that the 
chemical potential fi in the two (large ±u) unitary limits does 
not meet, nor does T c or A. This is because a (large) nega- 
tive u creates a filled impurity band below the main band, and 
is effective in bringing down the chemical potential, whereas 
a positive u creates an empty impurity band above the main 
band, and tends to raise /j, toward the particle-hole symmet- 
rical point, fi = 1. This result cannot and has not been ob- 
served in previous, non-self-consistent calculations where the 
real part of frequency renormalization is set to zero. 

In Fig. [T|, we compare the effect of the impurity concen- 
tration for different scattering strengths: the Born limit, both 
unitary limits (14 = ±100), and intermediate u = — 1. Both 
(b) T c and (a) A pg (T c ) are suppressed by increasing impurity 
density. This is natural in a model where the pseudogap origi- 
nates from incoherent pair excitations. As will be seen below, 
A P9 (T C ) is suppressed mainly because T c is lowered. Except 
in the Born limit, the chemical potential /i is fairly sensitive 
to rii, as shown in the inset. It is clear that the scattering in 
the Born limit is the most-effective in suppressing T c . In com- 
parison with experiment,El calculations at the AG level (i.e., 
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FIG. 13: Behavior of T c , fi (lower inset), and A (upper inset) as a 
function of —g/4t\\ in the unitary limit for u = —100 and rn = 
(solid curve), 0.02 (dotted), 0.05 (dashed), and 0. 1 (dot-dashed). For 
comparison, curves for 7 = 0.05 in the Born limit are also plotted 
(long-dashed). Impurity scattering in the Born limit is more effective 
in suppressing T c at relatively weak coupling. The parameters are 
n = 0.85, tx/t\\ = 0.01. 



the Born limit) tend to overestimate the T c suppression by as 
much as a factor of 2. This is in good agreement with the 
current result in the unitary limit. At large rij for large pos- 
itive it, the system is driven to the particle-hole symmetrical 
point, where the effective pair mass changes sign. It is usually 
hard to suppress T c by pairing at the particle-hole symmetri- 
cal point, as indicated by the solid curve in the lower panel. In 
fact, exactly at this point, the linear f2 term a' Q in the inverse T 
matrix expansion vanishes, so that one needs to go beyond the 
current approximation and expand up to the £1 2 term. Large 
negative u = —100 is more effective in suppressing T c and A 
than intermediate negative u = —1, in agreement with Fig. [IT] 
and the DOS shown in Fig. ||(a) and Fig. [k]. 

There is enough evidence that zinc impurities are attractive 
scatterer for electrons in the cupratesEJ Therefore, we con- 
centrate ourselves on negative u scattering in the unitary limit. 
Plotted in Fig. [j~3| are T c (main figure), /1 (lower inset), and A 
as a function of g for increasing n, = 0, 0.02, 0.05, 0.1 with 
u = — 100. Also plotted for comparison are the results assum- 
ing the Born limit with 7 = 0.05. Clearly, the Born limit is 
more effective in suppressing T c at relatively weak coupling, 
-ff/^n < 0.75, consistent with Fig. |l2| Both T c and A are 
suppressed continuously with m. However, it should be noted 
that a larger rii helps T c to survive a larger —5/4*11. This is 
mainly because the filled impurity band at large m pushes the 
system far away from particle-hole symmetry (see /i in the 
lower inset), and reduces the effective fermion density, so that 
the pair mobility is enhanced and the pair mass does not di- 
verge until a larger —5/4*11 is reached. 

To make contact with the cuprates, we use the non- 
double occupancy condition associated with the Mott insu- 
lator physics, as in Ref. 0, so that the effective hopping in- 
tegral is reduced to tu(x) & * x, where x = 1 — n is the 
hole concentration, and *q = 0.6 eV is the hopping integral 
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FIG. 14: Behavior of T c and A (inset) as a function of the hole dop- 
ing concentration x in the unitary limit for u — —100 and m = 
(solid curve), 0.02 (dotted), and 0.05 (dashed). The parameters are 
-g/Ata = 0.047, t±/t\\ = 0.003. 



in the absence of the on-site Coulomb repulsion. We assume 
—5/4*0 = 0.047, which is x independent. Then we can com- 
pute T c , [i, and A = A pg (T c ) as a function of x. The result for 
T c (main figure) and A (inset) is shown in Fig. [l4] for the clean 
system and m = 0.02 and 0.05. In the overdoped regime, T c , 
as well as the small A, are strongly suppressed by impurities. 
This provides a natural explanation for the experimental ob- 
servation that T c vanishes abruptly at large x; it is well-known 
that high crystallinity, clean samples are not available in the 
extreme overdoped regime. On the other hand, neither T c nor 
A is strongly suppressed in the highly underdoped regime, 
where the gap is too large. At this point, experimental data 
in this extreme underdoped regime are still not available. Our 
result about the suppression of T c and A in the less strongly 
underdoped regime ix-p- 0.1) are in agreement_with experi- 
mental observationsEH and other calculations. El 



It should be noted, however, that in our simple model, 
we do not consider the fact that disorder or impurities may 
reduce the dimensionality of the electron motion and thus 
suppress T c . Furthermore, since induced local spin and 
Kondo effects h&if been observed near zinc sites in both zinc=. 
doped YBCOBa and zinc-doped BisSrs-^La^CuOg+^El 
this raises an important question whether zinc can be treated 
as a nonmagnetic impurity. 



D. Gaps and superfluid density below T c in the presence of 
nonmagnetic impurities 

In this subsection, we study the effect of nonmagnetic im- 
purities on the behavior of the excitation gap A, the order pa- 
rameter A sc , and the pseudogap A pg as well as the chemical 
potential /iasa function of temperature in the superconduct- 
ing phase. The numerical solutions for these quantities are 
then used to study the temperature dependence of the super- 
fluid density n s /m at various impurity levels. We concen- 
trate on the unitary limit, which is regarded as relevant to the 



13 




FIG. 15: Behavior of the superfluid density n a /m and the excita- 
tion gap A (inset) in a d-wave BCS superconductor as a function 
of temperature T/T° for impurity concentration m = (clean, 
solid curve), 0.02 (dotted), and 0.05 (dashed) in the unitary limit 
u = —100. T° — 0.0416 is the T c in the clean system. The param- 
eters are n = 0.85, -g/At\\ = 0.3, and tx/t\\ = 0.01. 




FIG. 16: Behavior of the superfluid density n s /m and the various 
gaps (inset) in a d-wave pseudogapped superconductor as a function 
of temperature T/T® for impurity concentration m = (clean, solid 
curve), 0.02 (dotted), and 0.05 (dashed) in the unitary limit u = 
-100. HereT c ° = 0.0414. The parameters are n = 0.85, -g/^t\\ = 
0.5, andtxAy = 0.01. 



cuprates. To be specific, we use n = 0.85, t±/tn = 0.01, and 
u = —100 in the calculations presented below. 

We first study the impurity effect in the BCS case, without 
the complication of the pseudogap. Plotted in Fig. |l5| are the 
superfluid density n a /m (main figure) and the corresponding 
gap A in a <i-wave BCS superconductor as a function of the 
reduced temperature T/T c ° for the clean system (solid curves), 
impurity density m = 0.02 (dotted), and m — 0.05 (dashed) 
at -ff/4*|| = 0.3. Here T c ° = 0.0416 is the T c in the clean 
case. As expected, both T c and A(T), as well as n s /m, are 
suppressed by impurity scattering. In agreement with experi- 
ment, the low T normal fluid density is linear in T in the clean 
case, and becomes quadratic in the two dirty cases. 

Now we add pseudogap for the underdoped cuprates. We 
show in Fig. [16] the temperature dependence of n s /m (main 
figure) and various gaps (inset) in a d-wave pseudogapped su- 
perconductor for impurity concentration rii = (clean, solid 
curve), 0.02 (dotted), and 0.05 (dashed) in the unitary limit 
at —g/At\\ = 0.5. As the order parameter develops below 
T c , the pseudogap decreases with decreasing T. This reflects 
the fact that the pseudogap in the present model is a measure 
of the density of thermally excited pair excitations. The total 
gap A, the order parameter A sc , and the superfluid density 
n s /m are suppressed by increasing rii, similar to the BCS 
case above. However, at given T < T c , the pseudogap A. pg 
remains roughly unchanged. Furthermore, the low T power 
law for the superfluid density is different from the BCS case, 
as predicted in Eq. (|7|). It is now given by T + T 3 / 2 and 
T 2 + T 3 / 2 for the clean and dirty cases, respectively. Due 
to the presence of the T 3 / 2 term, the low T portion of the 
curves for n,; = 0.02 and 0.05 are clearly not as flat as in 
Fig. [fjl Nevertheless, it may be difficult to distinguish exper- 
imentally T 2 + T 3 / 2 from a pure T 2 power law. This T 3 / 2 
contributioruaf the pair excitations has been used to explain 
successfullyEI the quasi-universal behavior of the normalized 



superfluid densit*Ji s (T)/n s (0) as a function of T/T c . We 
have also founded preliminary experimental support for this 
T 3 / 2 term in low T penetration depth measurement in the 
cuprates as well as organic superconductors. Systematic ex- 
periments are needed to further verify this power law predic- 
tion. 

A careful look at the values of the zero temperature super- 
fluid density (n a /m)o f° r different values of rii in both Fig. [l5| 
and Fig. |l6| suggests that in the unitary limit, (n s /m)o drops 
faster with rii when rii is still small. This is manifested in a 
systematic study of (n a /m)o as a function of m, as shown 
in Fig. [n], with u = —100 (solid curve). This behavior has 
been observed experimentallyEJ Also plotted in the inset is the 
corresponding zero temperature gap Ao versus rii. Clearly, 
the slope d^/dn,i is much steeper as n, approaches zero, 
very different from the behavior of Ao. This demonstrates 
that n s /m is influenced more by the DOS than by the gap 
size. A very small amount of impurities may strongly sup- 
press (n s /m)o. This conclusion is significant in data analy- 
sis of the penetration depth measurement, especially when the. 
T = value of penetration depth is not measured directly.H 
For comparison, we also plot the corresponding curves in the 
Born limit. While the gap is suppressed faster, in contrast to 
unitary case, the slope d^f- / drii is smaller for smaller m. 



IV. DISCUSSION 

In Sec. II, we have used the approximation Eq. ( |23[ ) to bring 
the single-particle self-energy and thus the gap equation into 
a BCS-like form. This approximation derives from the diver- 
gence of the T matrix t pg (Q) as Q — > 0, which is the pairing 
instability condition. The spirit of this approximation is to 
"put" the incoherent, excited pairs into the condensate, by set- 
ting Q = 0. The contribution of these pseudo Cooper pairs 
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FIG. 17: Zero temperature superfluid density (n 3 /m)o and the gap 
Ao (inset) in a d-wave superconductor as a function of impurity con- 
centration rii in the unitary limit (it = —100, solid line) and in the 
Born limit (7 = m, dashed line). The parameters are n = 0.85, 
— p/4*n = 0.5, and t±/tu = 0.01. 



conductors. For s-wave superconductors, one would have to 
include self-consistently one more complex equation for the 
renormalization of Ak, when solvi ng fo r the renormalization 
of L). And the equations (37a) and (37c) also look much more 
complicated than Eq. ([39[). Nevertheless, since there is no 
node in the excitation gap for s-wave, the numerics is ex- 
pected to run faster. 

It is well-known that for d-waye superconductors, the An- 
derson's theorerro breaks downo For Anderson's theorem to 
hold, it requires that the gap and the frequency are renormal- 
ized in exactly the same fashion. This condition can be sat- 
isfied (approximately) only in weak coupling, isotropic BCS 
,s-wave superconductors, for which the real part of the fre- 
quency renormalization is negligible. Since the frequency u> 
is a scalar, this condition is violated when the gap Ak have 
any anisotropic dependence on k. Furthermore, when the gap 
is considerably large in comparison with the band width so 
that the upper limit of the energy integral cannot be extended 
to infinity, this condition will not be satisfied, either. In both 
cases, T c will be suppressed. 



to the single-particle excitation gap is calculated via Eqs. ( |l5| ) 
and (pl|), weighted by the Bose function. Therefore, the in- 
coherent pairs and the zero momentum condensate are not 
distinguished from each other in terms of the single particle 
self-energy, as they add up to a total excitation gap. However, 
they are distinct when phase sensitive quantities are involved, 
e.g., in the calculation of T c and of the superfluid density. 

With this approximation, there is a close analogy be- 
tween the Feynman diagrams in the current pairing fluctua- 
tion theory in Sec. II and its BCS counterpart in Appendix 
|a|. When the finite momentum pair propagators are removed 
(or "pushed into the condensate") from Fig. || through Fig. 
these diagrams will become their BCS counterpart in Fig. |l8| 
through Fig. (The diagram for BCS pairing vertex is not 
shown in Appendix |a|). 

This approximation is in general good when the gap is large 
in comparison with T, and when the impurity concentration 
is low. When the gap is small, the contribution of the inco- 
herent pair excitations is usually small, and does not have a 
strong effect on T c , so that T c is roughly determined by its 
BCS mean-field solution. When rii is large, the fermionic fre- 
quency renormalization is strong, and the pair dispersion also 
becomes highly damped. In this case, approximation Eq. ( p3f ) 
may not be quantitatively accurate. 

Even without the complication of impurities, incoherent 
pairs are not expected to deplete completely the spectral 
weight Wtbin the two quasiparticle peaks of the spectral 
function.Qo This, however, cannot be captured by the approx- 
imation Eq. (|23|). Unfortunately, we still do not know yet how 
to solve the Dyson's equations without this approximation due 
to technical difficulties. 

Another simplification comes from the <i-wave symmetry of 
the cuprate superconductors under study. Although we have 
kept the theoretical formalism general for both s- and <i-wave 
in Sec. II, the pairing vertex renormalization drops out when 
we finally carry out numerical calculations for d-wave super- 



V. CONCLUSIONS 

In this paper, we extend the pairing Auction theory to super- 
conductors in the presence of non-magnetic impurities. Both 
the pairing and impurity T-matrices are included and treated 
self-consistently. We obtain a set of three equations for (T c , fi, 
A(T C )) or (/i, A sc , A p£( ) at T < T c , with the complex equa- 
tions for the frequency renormalization. In consequence, we 
are able to study the impurity effects on T c , the order parame- 
ter, and the pseudogap. In particular, we carry out calculations 
for d-wave superconductors and apply to the cuprate super- 
conductors. Instead of studying the physical quantities with 
all possible combinations of the parameters rii, u, g, and n, 
we mainly concentrate on the negative u unitary limit, which, 
is regarded as relevant to the zinc impurities in the cuprates.EZl 

Calculations show that in addition to the low energy reso- 
nance in the imaginary part of the renormalized frequency, a 
considerably large | u | leads to a separate impurity band, with a 
spectral weight 2n,. The real part of the frequency renormal- 
ization, in general, cannot be set to zero in a self-consistent 
calculation. The chemical potential varies with the impurity 
concentration, so that the assumption of exact particle-hole 
symmetry is not justified when one studies the impurity ef- 
fects. One consequence of this chemical potential shift is that 
the repulsive and attractive unitary scattering limits do not 
meet as has been widely assumed in the non-self-consistent 
treatment in the literature. Unitary scatterers fill in the DOS 
mostly in the small uj region, whereas Born scatterers do in 
essentially the whole range within the gap. At small m and/or 
small u, there is a dip at uj = in the DOS, so that N(oj) van- 
ishes as a fractional power of u, which may in turn contribute 
a fractional power law for the low T temperature dependence 
of the penetration depth. 

Both T c and the pseudogap A pg (T c ) are suppressed by im- 
purities. In this respect, Born scatterers are about twice as ef- 
fective as unitary scatterer. Treating zinc impurities as unitary 
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FIG. 18: Impurity dressing at the Abrikosov-Gor'kov level in BCS 
theory. 



scatterers explains why the actual T c suppression is only half 
that predicted by calculations at the AG level (i.e., in the Born 
limit). In the overdoped regime, the gap is small, and there- 
fore the superconductivity can be easily destroyed by a small 
amount of impurities. In contrast, it takes a larger amount 
of impurities to destroy the large excitation gap in the under- 
doped regime. 

The reason Apg(T" c ) is suppressed is mainly because T c is 
suppressed. In fact, for a given T < T c , the pseudogap re- 
mains roughly unchanged (actually it increases slightly). The 
suppression of the total excitation gap arises from the suppres- 
sion of the order parameter. The density of incoherent pairs, 
as measured by n p , slightly increases for not-so-large m. This 
supports the notion that nonmagnetic impurities do not mainly 
break incoherent pairs—Instead, they scatter the Cooper pairs 
out of the condensate.LJ 

Our self-consistent calculations show that in the unitary 
limit, the low T superfluid density is quadratic in T in a BCS 
d-wave superconductor, in agreement with existing calcula- 
tions and experiment. Strong pair excitations add an addi- 
tional T 3 / 2 term, with preliminary experimental support. As 
a function of increasing rii, the zero T superfluid density de- 
creases faster at first for unitary scatterers, whereas the oppo- 
site holds for scattering in the Born limit. The former behavior 
is in agreement with experiment. 
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APPENDIX A: IMPURITY DRESSING FOR BCS THEORY 
AT THE ABRIKOSOV-GOR'KOV LEVEL 

In this appendix, we present the impurity dressing_for a 
BCS superconductor, following Abrikosov-Gor'kovESj, but 
in a more general form, namely, we do not assume G_ w = 
— G u . This will make it easier to understand the current the- 
ory in the presence of strong pairing correlations, as there is a 
strong similarity between the impurity dressing diagrams for 
both BCS theory and the pairing fluctuation theory. 



T 



FIG. 19: Replacement scheme from the AG level impurity scattering 
to self-consistent impurity treatment in BCS theory. 



For a pure BCS superconductor, we have the Gor'kov equa- 
tions, 

G -\K)G\K) = 1 - A k F ot (X) , (Ala) 



G° l (-K)F°\k) = AIG°{K) . 



(Alb) 



At the AG level, the relationship between the impurity 
dressed Green's functions G and F is represented by the Feyn- 
man diagrams shown in Fig. Define the impurity average 
G^j as in Eq. ©, and 



Ft^^Kk-k')| 2 Ff(J(') 



(A2) 



as well as their complex conjugate. Note Fig. 18 is actually 
Fig. 105 in Ref. [l4| Without giving details, we give the result 
following AG: 

( iuj _ £k _ G U )G(K) + (A k + F U )F\K) = 1 , (A3a) 



(iuj + e k + G-J)F\K) + (A k + Fl)G{K) = . (A3b) 

Define iuj = iuj — G u , iQ = —iuj — G- u , A k = A k + F u , 
and A k = A£ + Then we obtain 



G{K) 



F\K) 



iuj — e k 



(iQ - e k )(iw - ek) + A k A k 



K 

(iuj ~ e k )(iw - e k ) + A k A k 



(A4a) 



(A4b) 



For d-wave, the first equation becomes Eq. (^OJ). Note 
G(K) is no longer symmetrical in uj in general as a conse- 
quence of impurity scattering, but F(K) still is, since F(K) 
involves ±uo pairs. 

The above result can be easily extended to self-consistent 
impurity T-matrix calculations, by replacing the AG-level im- 
purity scattering with the self-consistent impurity T-matrices, 
as shown in Fig. [l^. The relationship between the regular 
and anomalous impurity T matrices T w and T^t are shown in 
Fig. One can easily write down the corresponding equa- 
tions, as follows. 



T u =u + uG^T^ - uFuT&\ 



(A5a) 
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FIG. 20: Relationship between impurity T matrices T u and T A t in 
BCS theory. 



T At = uFlTu + uG_ w T At 



Finally, one has 



u(l - uG-u) 



(A5b) 



(A6a) 



where G = £ k G ( K )> Ft = T,k F H K )> and similarly for 
their complex conjugat e. N ote t hese two equations are for- 
mally identical to Eqs. (37a) and (|37c|), except that the current 



T A t contains the factor A already. 



Now with the new definition iuj 
A k = A k + S A , and A* 
S w = n 4 T w and E A = n{T A u Eqs. ( |a1| ) for G and F remain 
valid. 



A* 



E A , as well as 
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